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ABSTRACT 

We propose that the majority of quasars at redshift z ~ 1 — 5 formed in the 
environment of new born collapsed halos with 1-D velocity dispersion ~ 400 km s -1 . 
The harboring coefficient / of quasars per halo and the lifetime of quasars depend 
only on local process, not modulated by the density inhomogeneities on scales larger 
than the size of the halos. Thus, the bias of quasars on scale larger than the size of 
these halos is mainly determined by the parameter a v used for quasar environment 
identification. With this model, the popular structure formation models, like SCDM 
and LCDM, can be fairly well reconciled with the data of quasars, including a. observed 
feature of the environment for quasars; b. redshift evolution of quasar abundance; 
c. the two-point correlation functions of quasars. This bias model predicts that the 
correlation function of quasars doesn't significantly evolve, or only slightly increases 
with redshift. 



Subject headings: cosmology: theory - quasars: general - large-scale structure of 
universe 



1. Introduction 



Mass distribution at high redshifts is being a hot subject of the large scale structure study. 
Data of various objects at moderate and high redshifts, in particular, clusters of galaxies, are 
becoming available for probing the formation and evolution of structures and for discriminating 
among popular dark matter models (e.g. Jing & Fang 1994; Eke, Cole & Frenk 1996; Bahcall, 
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Fan & Cen 1997; Kitayama & Suto 1997). Considering that quasars are the most distant among 
various luminous objects, they have also been applied in this approach (e.g. Bi & Fang 1997). 

However, as a mass tracer of the cosmic matter field, quasars are still playing a role different 
from clusters of galaxies. The problem stems from so-called "bias" . Clusters are a biased tracer 
of the mass distribution. The correlation amplitude of clusters is believed to be much higher 
than that of the underlying matter and increases strongly with the cluster richness (Bahcall & 
Soneira 1983). This bias is plausibly explained by the mechanism that the observed clusters are 
identified as massive collapsed halos of the density field (Kaiser 1984). That is, the bias of clusters 
is completely determined by the gravitational parameters, like mass M and virial radius r v ; Lr used 
to identify the halos. With this approach, a detailed confrontation can be made between theories 
and the observations of clusters. 

Quasars may also be biased tracer of the mass distribution. Recent observations indicate that 
the correlation amplitude of quasars may also be different from the underlying dark matter (Mo 
& Fang 1993; Romberg, Rravtsov & Lukash 1994; Croom & Shanks 1996; Franca, Andreani &; 
Cristiani, 1997). However, so far no detailed model is available for the bias of quasars, though 
their high clustering strength and environment imply that quasars are hosted by massive halos 
(see below). Because of the lack of such a model, one cannot confront the data of quasars with 
theoretical models in the way as for clusters. For instance, the abundance of quasars can only 
be used as an upper or lower limit to certain massive halos; no detailed comparison between the 
number densities of quasars and of halos is allowed. Obviously, it is very important to understand 
what kind of mass halos are associated with the majority of quasars. Such a knowledge will not 
only enable the observational data of quasars to be powerful tests for theoretical models of galaxy 
formation but also tell that what type of local environments is responsible for intriguing the 
nuclear activities of quasars. 

Like clusters and groups of galaxies, it is generally believed that quasars should be associated 
with certain type of collapsed dark matter halos. Yet, different from identification of clusters, the 
environment suitable for forming quasars is not merely determined by gravitational parameters, 
as the hydro processes are also involved. Therefore, the identification of quasar-harboring halos 
should be given by both gravitational and hydro parameters. In other words, not all halos with 
certain M and/or velocity dispersion a v harbor quasars, because certain hydro conditions must 
be satisfied. However, considering the hydro processes are local, it is reasonable to assume that 
the hydro conditions may not be modulated by the density inhomogeneities on scales much larger 
than the size of the halo I. In this case, the probability for a halo to have a quasar should be the 
same for all halos of the same kind, without depending on structures larger than I. Thus, the 
relative fractions of quasars with respect to the certain collapsed halos should be the same for all 
volumes larger than I 3 . Consequently, when averaged on scale larger than I, the distribution of 
quasars n qso (r, z) at redshift z should be proportional to that of the considered halos, n^ (r, z). 
Thus, all effects of the hydro processes can be absorbed into a "normalization factor" A, i.e. 
n qso (r, z) = Arih a i (r, z), and A is less dependent on z than n^ D (r, z). The bias of quasar 
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distribution with respect to the mass distribution is then dominated by the bias of the selected 
halos with respect to the mass. Based on this analysis, quasar bias, at least on large scales, 
may also be only gravitational, depending on the gravitational parameters used for selecting the 
quasar-suitable halos. 

Accordingly, a possible model for quasar bias should at least satisfy the three conditions. 1. 
Gravitational environment given by the identified halos is consistent with the observed environment 
around quasars; 2. The abundance of quasars, n qso (z), at redshift z is proportional to the number 
density of the identified halos, nh a io( z ) m a redshift-independent way, i.e. n qso {z) = An^aio where 
A is a z-independent constant, 3. The amplitude and z-evolution of the halo-halo correlation 
function are consistent with the observed correlation function of quasars. In this letter, we will 
show within the framework of the CDM cosmogonic theories that such a bias model can indeed be 
settled following the above-mentioned points. The details of the points 1, 2 and 3 will be discussed 
in the §2, 3 and 4, respectively. 

2. Basic assumption: The gravitational condition for a quasar halo 

A bias model of quasars is just a phenomenological relationship between the cosmic density 
field and quasars. With this relationship quasars can be identified from the mass density field. In 
this sense, bias, in fact, is a model for the environment suitable for the quasar formation. 

In what environment are quasars most likely to be formed? Many observations indicate that 
quasars are preferentially located in groups of galaxies. The evidences include the quasar-galaxy 
covariance function (Yee & Green 1987), the galaxy environments around quasars (Ellingson, et 
al. 1991), clustering of quasars (Bahcall & Chokshi 1991; Mo & Fang 1993; Komberg, Kravtsov & 
Lukash 1994; Croom & Shanks 1996; Franca, Andreani & Cristiani, 1997), and the CIV-associated 
absorption in high redshift radio-loud quasars (Foltz et al 1988). Recently, an observation of 
a companion to quasar BR1202-0725 with high redshift z = 4.7 directly shows that the width 
of their Lya emission lines is around 400 km s _1 (Petitjean et al. 1996.). These observations 
seemingly point to quasars being identified with the newly collapsed halos with 1-D velocity 
dispersion like groups, i.e. a]f ~ 400 km s -1 . It should be pointed out that this environment 
condition may not be necessary for low redshift quasars (Smith, Boyle & Maddox, 1995), because 
galaxies and clusters underwent a significant evolution at z ~ 0.5. But, it is reasonable to assume 
that the environment with o\ d « 400 km s -1 is favored by the formation of quasars at higher 
redshifts: a) enough collapsed objects to form the engine of a quasar; b) dispersed gas to feed 
the engine through accretion; and c) not too many proto-galaxies to possibly disrupt the process 
of quasar formation. In the next two sections, we shall compare, within the framework of CDM 
cosmogonic models, the theoretical predictions for such halos with the observed quasar abundance 
and correlation functions. 
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3. Redshift Evolution of Abundance of quasars 

For a Gaussian initial perturbation, the comoving number density of halos with 1-D velocity 
dispersion a\ d can be calculated with the Press & Schechter formalism (1974) as 

d\nA(R,z) S c ( S 2 C 




" dlnR A(R,z) eW \ ~A 2 (R,z 

where R is the Lagrangian radius of the dark halo being considered and S c ~ 1.69 almost 
independent of cosmologies. A(R, z) is the rms of the linear density fluctuations at redshift z 
within a top-hat window of radius R, and is determined by the initial density spectrum P(k) 
and normalization factor a§ = A(8 /i _1 Mpc, 0). The relationship between a]f and R is given by 
(Narayan k White 1988) 

al d = c^ 2 HMl + z) 1 ' 2 . (2) 

for a universe with the density parameter Qq. The N-body simulation results (Jing & Fang 1994) 
showed c a ~ 1.2 which is used for all calculations in this work. Two representative CDM models, 
i.e. the standard CDM (SCDM) and flat low-density CDM (LCDM) models are employed. The 
Hubble constant h = i^o/lOOkms" 1 Mpc" 1 , mass density £lo, cosmological constant Ao and as 
are taken to be (0.5, 1,0,0.58) and (0.75,0.3,0.7, 1) for the SCDM and LCDM, respectively. The 
models with these parameters seem to provide a good approximation to many observational 
properties of the Universe, especially the abundance of clusters (which is much related to the topic 
of this work). 

The total number density of the collapsed halos with the velocity dispersion greater than a 
certain value, say o\f m -, is 

N(>alf m ,z)= n(a 1 v d ,z)da 1 v d . (3) 
iff" 

The birth rate of halos with a]f > a\ d m is dN(> aj d n ,z)/dt. This birth rate is shown in Fig.l 
for a\ d m = 200, 400, and 800 kms -1 . For each a} d m , the birth rate in the two models possesses 
similar shape. At <r}f m = 400kms _1 , the birth rate keeps steady from z = 5 to about 2, and then 
rapidly drops to zero at z ~ 0.3 for SCDM and 0.7 for LCDM. The peak birth rate however shifts 
to higher redshift for smaller a\f m , and to lower redshift for larger a\ d m . As has been discussed 
in §1 and 2, each newly collapsed halo with a\ d > a] d m may host / quasars in average, and the 
harboring coefficient / is ^-independent. If the mean lifetime of quasars is t qso , which is also 
z-independent, the comoving number density of quasars at the epoch of redshift z is given by 



, , , /■*« dN(>a} d n ,z')dz^ 
n{z) = f / - l m ' —dt. 4 

Jt(z)-t aso dzf dt 



An accurate quasar lifetime t qso is not important only if it is z-independent and is much less than 
the age of the universe for redshift z < 5. Considering that birth rate of the halos is slowly varying 
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with redshift when z > 1, we have approximately n(z) ~ ft qso [dN(> (r}f m , z)/dz](dz/dt). Since 
the shape of the birth rate is a strong function of the velocity dispersion (especially as a function 
of z), the abundance of quasars provides a strong test to the bias model proposed here. 

Fig.2 plots n{z) for the SCDM and LCDM models. The "normalization constant" t qso f are 
adjusted in order that the theoretical maximum number densities of quasar abundance can fit 
with the observed one. The data points are the number density of quasars complete to absolute 
magnitude M B = -26 (Pei, 1995, Hewett, Foltz & Chaffee 1993, Schmidt, Schneider k Gunn 
1992). In these observations, the density n(z) is measured in the Einstein - de Sitter cosmology 
(f2 = 1, Ao = and h = 0.5). These data have been corrected to the case of non-zero Ao in 
the panel of LCDM. The figure shows the redshift evolution of quasar abundances is fairly well 
described by the newly collapsed halos with ajf m = 400 km s -1 . Both the SCDM and LCDM 
models are in reasonably good agreement with the observed abundance n qso (z) if the constant 
(z-independent) parameter t qso f is taken to be 0.02- 10 7 yr for the SCDM model and 0.33- 10 7 yr for 
the LCDM model, though the best fitting value of alf is slightly lower (about 340 km s -1 ) for the 
SCDM model. However for a\f m = 200, and 800 kms -1 , the predictions cannot fit with observed 
redshift evolution regardless of how to adjust the "normalization constant" t qso f. Therefore, 
the consistency between theoretical and observed redshift evolution of quasar abundance can be 
achieved only for <j]f ~ 400 kms -1 halos, which is not trivial. This result strongly supports that 
the majority of quasars be associated with new born halos with 1-D a\f m ~ 400 kms -1 . 

Particularly, for an environment with a given mass or velocity dispersion, the luminosities 
and lifetimes of quasars are still dispersed. However, if the distributions of the luminosities and of 
the lifetimes for this type of halos are not modulated by large scale perturbations, all conclusion 
here should hold. 



4. Amplitude and redshift-evolution of quasar correlation function 

According to the analysis of §1, the bias model for a\f m ~ 400 km s -1 halos is available on 
the scales larger than their typical size which is about 3(1 + z)~ 1 / 2 Q^^li' 1 Mpc. Therefore, 
on scales larger than 5 h" 1 Mpc the bias of quasar distribution is due mainly to the bias of a]f 
halos. On such large scales and at high redshifts, the mass distribution is still linear. Therefore, 
the two-point correlation function of the halos at redshift z is given by 



£(r, z) =b 2 (R,z)U(r,z) 



(5) 



where r is physical radius, £ m (r, z) is the linear mass correlation function, and the bias factor b is 
given by (Mo & White 1996) 



b(R,z) = 1 + 



A 2 (R,z) 



1 



(6) 
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Fig.3a shows the correlation functions for a\^ m = 200, 400 and 800 kms halos at z = 2. It 
shows that for o\f m = 400 km s -1 the correlation length ro defined by £(r) = (r/ro)~ L8 is ~ 6 h _1 
Mpc for the SCDM, and ~ 14 h _1 Mpc for the LCDM. Observationally, the two-point correlation 
function of quasars is found to obey the same power law. With qo taken to be 0.5, the correlation 
length r at z = 1.5 is found to be ~ 6.6 ± 0.5 h' 1 Mpc (Mo & Fang 1993), 10 ± 2 h~ l Mpc 
(Komberg, Kravtsov & Lukash 1994), and 6.2 — 8.0 hr 1 Mpc (Franca, Andreani & Cristiani 
1997). Actually, these results do not sensitively depend on qo, because on the length scales of 
about 10 h~ l Mpc, the influence of the parameter q$ is small. We should be careful in comparing 
the observed data with Fig. 3a. Generally, the samples of quasars employed for the correlation 
statistics possess different limit magnitude for different redshifts. The samples employed for the 
correlation statistics contain quasars with Mj, > —26. Namely, they are not complete in the sense 
as the data used in Fig. 2. The larger Mb may lead to a lower amplitude of the correlations. 
Considering this uncertainty, both SCDM and LCDM with o\f m = 400 km s -1 halos are acceptable, 
though a smaller or larger value of o\£ m may also be tolerated by the clustering observation alone. 

The redshift evolution of the halo-halo correlations for SCDM and LCDM is plotted in Fig.3b. 
It is interesting to see that the amplitudes of the correlation functions do not significantly evolve 
with redshift, having only a slight increase with redshift. This is because the clustering in the 
mass distribution of dark matter £ m (r,z) always increases with time, or decreases with redshift, 
but the bias factor b(R, z) increases with redshift. The two effects seem to be balanced by each 
other, giving a very slowly varying of £(r, z) with z. 

The z-evolution of quasar clustering has been studied for more than one decade (e.g. Chu & 
Fang 1987, Shaver 1988) and the results are quite scatter. Some showed a weak decrease of the 
correlation amplitudes around redshift 1.5 on scales of > 10 /i _1 Mpc (Mo & Fang 1993; Komberg, 
Kravtsov & Lukash 1994; Croom & Shanks 1996.) Some concluded no significant z-dependence 
from z ~ 1.5 — 2.9 (Zitelli et al. 1992.) And some even show very weak increase with redshift. 
Recently, a weak z-increasing correlation from z < 1.4 to > 1.4 is reported (Franca, Andreani & 
Cristiani, 1997.) Obviously, these diverse data cannot provide a concrete test on the prediction 
of Fig.3b. Although there is quite a bit uncertainty in the observed redshift evolution, the 
current result - no significant evolution of either z-increase and z-decrease - is consistent with the 
developed bias models. 



5. Conclusions 

We showed that velocity-dispersion-selected halos are a possible mechanism for the bias of 
quasars. The majority of quasars at redshift z ~ 1 — 5 formed in the environment of new born 
collapsed halos with 1-D velocity dispersion a]f" ~ 400 kms -1 . Both the harboring coefficient 
/ per halo and the lifetime of quasars are z-independent. With this bias model, the popular 
structure formation models, like SCDM and LCDM, can be fairly well reconciled with data of the 
abundance and correlations of quasars at z > 0.5. 
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It is interesting to point out that the velocity dispersion identified halos generally don't have 
the same mass. Eq.(2) shows that for a given cr^f, the redshifts the higher, the mass of the halos 
the smaller. This result has already been recognized in an earlier study, which shows that in order 
to fit with quasar abundance at high redshift, the mass of the halos has to be smaller than at the 
lower redshift (Bi & Fang 1997). 

With this model, one can predict that 1. The environment for quasars at redshifts from z ~ 1 
to 5 should be characterized by a velocity dispersion, a^ d ~ 400 km s -1 ; 2. The amplitudes of 
quasar two-point correlation function at high redshifts don't significantly evolve with redshifts. In 
the paper, only the models of the SCDM and LCDM are considered. We can expected that with 
better data of quasars becoming available, the bias model of quasars will play more important role 
for discriminating among models of structure formations. 
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Fig. 1. — The birth rate of halos with the 1-D velocity dispersion larger than o\f m . From top to 
bottom curves, a\f m are 200 km s -1 , 400 km s -1 , and 800 km s -1 respectively. 
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Fig. 2. — The quasar mean density predicted by the model (solid curves) versus the observation 
data (dots; see text). The curves from right to left are for (j\f m = 200 km s^ 1 , 400 km s -1 , and 
800 km s -1 respectively. The mean lifetime t qso for quasar is taken to be 10 7 yr. The harboring rate 
/ is adjusted to match the theoretical maximum density to the observed one. In the increasing 
order of aff m , f is 0.45 x 1(T 3 , 0.19 x 10" 1 and 2.1 for the SCDM, and 0.84 x 10" 2 , 0.33, 38.6 for 
the LCDM. 
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Fig. 3. — a) The correlation function at z = 2 of the halos with the 1-D velocity dispersion equal to 
(j\f m . The curves from bottom to top are for a}f m = 200 km s -1 , 400 km s -1 , and 800 km s -1 ; — b) 
The evolution of the halo correlation function on the linear scale. The halo correlation function is 
normalized by the mass correlation function at z = 0. The 1-D velocity dispersion for the halos is 
taken to be 200 km s" 1 , 400 km s" 1 , and 800 km s -1 (curves from bottom to top). 



